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We study the effects of surfaces on resistor percolation at the instance of a semi-infinite geometry. 
' Particularly we are interested in the average resistance between two connected ports located on 

the surface. Based on general grounds as symmetries and relevance we introduce a field theoretic 
Hamiltonian for semi-infinite random resistor networks. We show that the surface contributes to the 
£f} ' average resistance only in terms of corrections to scaling. These corrections are governed by surface 

resistance exponents. We carry out renormalization group improved perturbation calculations for 
the special and the ordinary transition. We calculate the surface resistance exponents (f>s and ^g 3 
for the special and the ordinary transition, respectively, to one-loop order. 
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to ' I. INTRODUCTION 



Percolation |jj is the perhaps simplest model for the irregular geometry which occurs in disordered media. Per- 
colation is intuitively appealing and it has a large variety of applications. Moreover, percolation is the prototype 
of a geometric phase transition. The continuing interest in percolation results in an abundance of publications year 
c~j . after year. In particular the purely geometric aspects of percolation have been studied extensively. From the current 
perspective, non equilibrium properties of percolation, like transport on percolation clusters, are of growing interest. 
I O, , Of course these issues are typically more challenging than the equilibrium properties. 

Random resistor networks (RRN) play a major role in the study of transport on percolation clusters. A RRN is 
fSJ , simply a bond percolation model in which the occupied bonds are assigned a finite, nonzero conductivity. Commonly 
t> ■ studied in the theory of RRN is the average resistance Mr(x,x') between two connected ports x and x' when an 
external current is inserted at x and withdrawn at x'. It was found that Mr(x, x') scales at the percolation point 
. as 

£j : M R {x,x')~\x-x'f /v , (1.1) 

with (j> being referred to as resistance exponent and where v is the correlation length exponent of the percolation 
, universality class. 

The renormalization group provides a powerful and elaborate framework to investigate RRN analytically. In 
particular a field theoretic approach based on the seminal work of Stephen Q and Harris and Lubensky j|] has 
1 ' prooved to be fruitful. Using this approach the resistance exponent has been calculated || to second order in 
e = 6 — d, where d is the spatial dimension. Moreover, the approach has been used to compute <f> for continuum 
percolating networks |7|, cf> for diluted networks of nonlinear resistors [js] -|To|, several fractal dimensions characterizing 
percolation clusters |p-icf , an entire family of multifractal exponents for the moments of the current distribution in 
RRN etc. 

Though a field theory of boundary critical phenomena has been established (for background on the field theoretic 
approach to boundary critical behaviour see Refs. 0,^5)) and successfully applied to geometric percolation |L6 17 
j_J the field theoretic approach to RRN has not yet been extended to include surfaces. In this paper we present such an 
extension based on the approach of Stephen and Harris and Lubensky. We consider the, at least from the standpoint of 
field theory, simplest geometry which comprises a surface, viz. a semi-infintc geometry. The central question addressed 
by this paper is: What is the critical behavior of Mr(x,x') when the ports x and x' are located on the surface and 
how does the surface contribute to it? 
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The plan of presentation is the following: the paper has two main sections, Sees. [n] and HI. In Sec. [n] we provide 
some background on the phenomenology of resistor percolation. Then we develop a field theoretic model which has its 
manifestation in a Landau-Ginzburg- Wilson type Hamiltonian. We show that the term in this Hamiltonian steaming 



we 



from the conductance of the surface bonds is irrelevant in the sense of the renormalization group. In Sec. II I 
present the core of our renormalization group analysis. We calculate correction to scaling exponents for the average 
resistance which are associated with the irrelevant surface term. These are then compared to the Wegner exponent 
for percolation. Section |y| contains our conclusions. Details of the calculations are relegated to three appendices. 

II. THE MODEL 
A. Semi-infinite random resistor networks 

Consider bond percolation on a semi-infinite lattice in d dimensions bounded by a [d — l)-dimensional plane. A 
lattice site i is either located off the surface, i.e., i € B = {x = (x\\,z) \ x\\ & Z d ~ x ,z £ N/{0}} , or on the surface, 
i.e., i <E S = {x = (x\\, 0) | x\\ € Z d_1 } . Each bond between nearest neighbors on the surface is occupied by a resistor 
of conductance ag with probability p$ or unoccupied with probability 1 — p$. The bonds between all other nearest 
neighbor pairs are occupied by resistors of conductance a with probability p or empty with probability 1 — p. 

Different phases can be distinguished depending on the values of p and ps [0 • A sketch of the phase diagram 
is given in Fig. [l]. The parameter r ~ p c — p indicates if p exceeds its critical value p c . ps may be enhanced with 
respect to p. The parameter c ~ ps, c (p = Pc) — Ps serves as a measure of the enhancement. Let P per c (x) denote the 
percolation probability that site x belongs to an infinite cluster. Assume c > 0. Crossing over from r > to t < 
leads to the formation of an infinite cluster. However, P porc (x € S) < P pcrc (x <E B) by virtue of the missing neighbors 
at the surface. The phase transition that takes place at r = is called ordinary transition. Now suppose that r is 
subcritical and that the system crosses over from p$ < ps,c to ps > Ps,c at fixed p. At this so-called surface transition 
an infinite cluster forms in the vicinity of the surface but P pe rc (( x \\,z)) fa^s off exponentially for z — » oo. Upon 
increasing p at fixed ps, the surface transition is followed by the extraordinary transition at which the exponential 
decay ceases to exist. The point (c = 0, t = 0) defines a tricritical point describing the so-called special transition. 

Suppose a current / is injected into a cluster at site x and withdrawn at site x' . The current carrying bonds 
constitute, apart from Wheatstone bridge type configurations, the backbone between x and x' . The power dissipated 
on the backbone is by definition 

P = I(V X -V X >) , (2.1) 
where V x is the potential at site x. Using Ohm's law, 

(Vj - Vi) = I id , (2.2) 
where p j is the current flowing through the bond from j to i, it may be expressed entirely in terms of voltages as 

P = R&x')- 1 (V x ~ V x ,f = Y, *U (V t - V 3 f = P({V}) . (2.3) 

{id) 

Here R(x,x') is the total resistance of the backbone, the sum is taken over all nearest neighbor pairs on the cluster 
and {V} denotes the corresponding set of voltages. As a consequence of the variation principle 
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^({VB-E 1 ^ 



= , (2.4) 



one obtains Kirchhoff 's law 



Y J <TiAVi-V j )=-Y,Iij = Ii , (2.5) 

U) (j) 



where h — I (5i, x — Sj, x ' ) an d the summations extend over the nearest neighbors of i. 
Alternatively to Eq. ([2.3]) the power can by rewritten in terms of the currents as 



P = R(x, x')I 2 = J2 Pi^h = P (W) ' ( 2 - 6 ) 

{id) 
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with {/} denoting the set of currents flowing through the individual bonds and pij — a { j. Obviously the cluster may 
contain closed loops as subnetworks. Suppose there are currents {-f^} circulating independently around a complete 
set of independent closed loops. Then the power is not only a function of / but also of the set of loop currents. The 
potential drop around closed loops is zero. This gives rise to the variation principle 



d 



l) =0. (2.7) 



Eq. (2.7) may be used to eliminate the loop currents and thus provides us with a method to determine the total 
resistance of the backbone via Eq. ( |2.6| ) . 

Since the resistance of the backbone depends on the configurations C of the randomly occupied bonds, one introduces 
an average (• • -)c over these configurations. It is important to recognize that the resistance between disconnected sites 
is infinite. Therefore one considers only those sites x and x' known to be on the same cluster. Practically this is done 
by introducing the indicator function x(x,x') which, for a given configuration C, is unity if x and x' are connected 
and zero otherwise. Then the nth moment of the resistance R with respect to the average (• • -)c subject to x and x' 
being on the same cluster is given by 

m£> = ( X (x, x')R{x, x'T) c I (X(x, x')) c . (2.8) 



B. Generating function 



Our aim is to determine the average resistance Mr — M„ . Hence our task is twofold: we need to solve the 
set of Kirchhoff 's equations ( |2.5| ) and to perform the average over all configurations of the diluted lattice. It can be 

accomplished by employing the replica technique . The voltages are replicated Z?-fold: V x — ► V x = (Vx^ , • ■ ■ , vj 1 ^ 

One introduces 



i>%(x) = exp (iX- V x ^j , 
where A • V x = J2 a A < - Q - ) V^ a ' ) and A ^ 0. The corresponding correlation functions 

G(x,x';\)=(iP x (x)iP_ x (x')) iep 

are defined as 

G (x, A) = hm /z-° J UU dV ? ^ {~\ P ({?}) + f E ^ + ^ ■ {?. 

\ j Q=l \ I 

Here P ({^}) = Ei,j, a <*i,j {yl^ - ^/ a) ) and z is the normalization 




(2.9) 



(2-10) 



(2.11) 



Z = J]JdVie X p 



(2.12) 



Note that we have introduced an additional power term ^ J2i V? ■ This is necessary to give the integrals in Eqs. (2.11) 



and ( 2.12 ) a well defined meaning. Without this term the integrands depend only on voltage differences and the 
integrals are divergent. Physically the new term corresponds to grounding each lattice site by a capacitor of unit 
capacity. The original situation may be restored by taking the limit of vanishing frequency, lo — > 0. 

The integrations in Eq. ( [2.1l| ) can be carried out by employing the saddle point method. Since the integrations are 
Gaussian the saddle point method is exact in this case. The saddle point equation is identical to the variation principle 
stated in Eq. (2.4). Thus the maximum of the integrand is determined by the solution of Kirchhoff's equations ( |2.5| ) 
and 



G ( x, x': A 



exp i-—R(x,x') 



(2.13) 
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up to an unimpo rtant multiplicative constant which goes to one in the limit D — > 0. Taylor expansion of the right 
hand side of Eq. (|2.13|) about A (q) = leads to 



G 



A 2 



(x, x'; X) = (x(x, x')) c I I - Y M R ] ( x > x ') + 



1 

k\ 



A" \ . ( fc) 



(as, as 7 ) + 



(2.14) 



Hence the correlation function G can be exploited as a generating function for the moments Md given by 

<9" 



( X (z, z'))c m£° (a:, x 1 ) = - G (s, a/; X) 



a(-A 2 /2) 



A=0 



(2.15) 



C. Field theoretic Hamiltonian 



Since infinite voltage drops between different clusters may occur, it is not guaranteed that Z stays finite, i.e., the 
limit limjj^o Z D is not well defined. Moreover, A = has to be excluded properly. Both problems can be handled by 



resorting to a lattice regularization of the integrals in Eqs. (2.11) and (2.12). One switches to voltage variables 9 = A9k 
taking discrete values on a D-dimensional torus, i.e. k is chosen to be an £>-dimensional integer with — M < k^ < M 
and k^ = fcWmod(2M). A9 = 9m /M is the gap between successive voltages and 9m is the voltage cutoff. In this 
discrete picture there are (2M) D — 1 independent state variables per lattice site and one can introduce the Potts spins 



- (2M)- D £ exp (*A ■ 0) ^ x (x) = 5^ - {2M)-° 



(2.16) 



A#0 



subject to the conditio n J2a ®a( x ) ~ 0- 

Now we revisit Eq. (2.11). Carrying out the average over the diluted lattice configurations provides us with the 
weight exp(— H icp ) of the average (...) rep : 



, / 1 iuj 

H Icp = - In ( exp I --P + — 



E ln ( exp ( ~2° id ~ ® j ) 



E 



IU1 

~2 



E 



a2 



(2.17) 



By dropping a constant term Ng ln(l — p) + N$ ln(l — ps) with N$ being the number of bonds between sites on the 
surface of the undiluted lattice and TVg being the corresponding quantity related to remaining bonds we obtain 



i£B,jeBl)S 



ies,jes 



= E EM^W*) E E^^-^WCO^C?) 



ieBjeBus i 



ies,jes , 



EEK*K (i) ■ 



(2.18) 



where 



and 







= lnO 



1-J» 



exp ——<j9 



(2.19) 



^ 9 =ln 1 



1 - ps 



exp [ — • ^ 



(2.20) 
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The summations in Eq. ( 2.1E ) run over nearest neighbor pairs subject to the specified conditions. Note that Kq {d^ 

and Kg (Oj are exponentially decreasing functions in replica space with a decay rate proportional to er -1 and a^ 1 
respectively. The Hamiltonian H Tep describes interaction of Potts spins in a semi-infinite system with an external one 
site potential h (jij . The interaction is rotationally and translationally invariant in replica space. For large a and as 

the interaction is short ranged not only in real but also in replica space. Moreover, the interaction potential K B 
is an analytic function of 9 2 . Thus the Fourier transform 



K B (X) 



J2 ex P In ( 1 + -!— exp ( --a6 ■ 6 \ ) (2.21) 



(2M) D ^""^ \ "" "J~\ • 2 

9 

can be Taylor expanded as 

oo 

K B (X) =r + ^ Wp (X 2 ) P , (2.22) 
p=i 

with t and w p ~ a~ p being expansion coefficients. Analogously one obtains for the Fourier transform of K$'- 

oo 

K s (X) = c + J2 w s., P (X 2 ) P , (2.23) 

P =i 

with coefficients c and ws, P ~ c 5 p . 

In the limit of perfect transport, a — > oo and erg — ► oo, i^g and iCg go to their local limits Kg (o^j = KgSg^ 

and Kg (o^j = K$5gQ, with Kg and Kg being positive constants. The interaction part of the Hamiltonian reduces to 
H% = -K B J2 E ®e W *g U)-Ks £ ]T *g (i) $ fl - (i) . (2.24) 

leBJeBUS g i&B,j&B g 

This represents nothing more than the semi-infinite (2M) D states Potts model which is invariant against all (2M) D \ 
permutations of the Potts spins <f>g. If a^ 1 =/= and a^ 1 ^ 0, this S^ 2 m) d symmetry is lost in favor of the short range 
interactions. 



We proceed with the usual coarse graining step and replace the Potts spins <&g(x) by order parameter fields <p ^x. 

which inherit the constraint ^2g<fi ^x, (fj = 0. We model the corresponding field theoretic Hamiltonian H in the spirit 

of Landau as a mesoscopic free energy from local monomials of the order parameter field and its gradients in real 
and replica space. The gradient expansion is justified since the interaction is short ranged in both spaces. Purely 
local terms in replica space have to respect the full S^ 2 m) d Potts symmetry. After these remarks we write down the 
Landau-Ginzburg- Wilson type Hamiltonian 

H = H B +H S , (2.25a) 

with the bulk contribution 

n B = J v d d x E ( x ' Q) K ® ( A < A e) V ( x ' 0) + | <P ( x > ^) 3 } . ( 2 - 25b ) 



and the surface part 

Us = \^Ps (x,| , 9) K S (Ag) ^ s (x,| , 6) \ , (2.25c) 

6 

where 

K B (A, Ag) = t - A - wAg (2.25d) 
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and 



(2.25e) 



The integration in Eq. ( 2.251:) extends over the half-space V — {x = (x|| , z) | X|| e R d_1 ,z > 0} whereas the 
integration in Eq. (2.25c ) ex tends over the corresponding boundary dV . Terms of higher order in the fields have 
been neglected in Eq. ( 2.25b| ) and Eq. (2.25c) since they turn out to be irrelevant in the rcn ormalization group sense. 
Moreover, a term of the form ipsd n (p = lim 2 ^ + ¥><s3 z </? has been neglected in Eq. (2.25c) though it is marginal on 
dimensional grounds. Such a term turns out to be redundant, cf. Ref. p^ j. r, c, w, and ws are now coarse grained 
analogues of the original coefficients. Note that, u pon setting w = ws =0,% reduces to the usual semi-infinite Potts 
model Hamiltonian as studied by Diehl and Lam n% . 



D. Irrelevance of ws 

Here we discuss the relevance of ws in the renormalization group sense. Let P denote the set of parameters 
{t, w, c, ws} and b some scaling factor for the voltage variable: 9 — > bO. By substitution of ip ^x, e) = ip' ^x, b9) and 

ips ( x, 9 ) = ip' s ( X|| , b9 ) into the Hamiltonian we get 



H 



<f' x,66» ,<p' s I X|| 



P 



= Jd d xJ2{ \p' (x, bO) K B (A, A e -) <p? (x, , 

V 8 

+ J dv dd ~ lx W E { \v's (x|| , bO) K S (A,-) & (x,| , W 



(2.26) 



Renaming the scaled voltage variables, 9' = b9, leads to 



n 



f |X,C/ ,<f s I X| 



9') ,P 



J d d x^ \\<p' (V) K b (A,6 2 A e -,) V ' (x,0') + |y/ (x,^ 



(2.27) 



Clearly a scaling of the voltage variable results in a scaling of the voltage cutoff: 6m — > However, by taking the 

limit D — > before 0m — * oo the dependence of the theory on the cutoff drops out. We can identify 9' and 9 and thus 



H 



n 



<pU,e) ,<p s U h e) ,p' 



ip fx, be) , (p s (xn ,be) ,p 

where P' — {r,b 2 w,c,b 2 ws} ■ Equation ( p. 28 ) implies for the connected correlation functions 

IN m \ conn 

G n ,m N x, X|| , e\ ;r, w, c, w 5 ) = hm ( J| V ( x i, J| <£s ( x |b , ^j) > 

^° \»=i j=i V / w 

that 

Gjv.M f{x, X|| , e \ ] t, w, c, w s *) = G n .m ( |x, X|| , b0 j ; t, 6 2 w, c, 6 2 



(2.28) 



(2.29) 



(2.30) 



Note that x in Gn,m refers to a point located off the surface, i.e, x = (x||, z) with z > 0. In other quantities, however, 
x may refer to an arbitrary point. The two-point correlation functions Gn,m with N + M = 2 are nothing more than 



replica space Fourier transforms of ( ■0r(x)-0_r(x) ) . We deduce from Eq. (2.14) that 



n 



G 



X 2 M R ((x, x') ; r, w, c, w s ) = (b^Xj M R ((x, x') ; r, b 2 w, c, 6 2 u; 5 ) . (2.31) 

The freedom to choose b has not been exploited yet. With the choice b 2 = w~ x the previous scaling relation turns 
into 

M R ((x, x') ;t,,w,c,w s ) = wf ((x, x') ; r, c, ^) , (2.32) 

where / is a scaling function. Dimensional analysis of the Hamiltonian shows that w\ 2 ~ /i 2 and ~ where 

H is the usual inverse length scale. Thus ws/w ~ /i , i.e., ws is an irrelevant coupling. 



III. RENORMALIZATION GROUP ANALYSIS 



The three types of transition occurring at bulk criticality r = 0, namely the ordinary, the special, and the ex- 
traordinary transition, are described by renormalization group (RG) fixed points where c takes the value c* rd = 00, 
c* p = 0, and c* x = —00 respectively (cf. Sec. IIIC3 of Rev. O] and references therein). In the remainder of this 
paper we focus on the special and the ordinary transition. This section provides an outline of our renormalization 
group improved perturbation calculation. 



A. Gaussian propagator 



Now we determine the Gaussian propagator for our model. For this purpose the irrelevant term has to be discarded 
from the Hamiltonian. Then the saddle point solution of 

J V<pcxp(-H) , (3.1) 
for all x and 8, is determined by the so-called 



where Dip indicates an integration over the set of variables |</? ^x, | 
equations of motion 

Tip — Aip — wAgip + ^-ip 2 = 

and 



(3.2a) 



Ps+d n p s = Q. (3.2b) 

In Eq. ( 3.2b ) d n ps ( x |h^) denotes the normal derivative lim z _> + d z tp ^x, 9J . As a consequence of Eq. ( |3.2| ), the 
Fourier transform G p % (z, z') of the Gaussian propagator G (xy — xji, z, z' , — 6'j has to satisfy 



t + p 2 + wX 2 - d 2 



G prx {z,z')=5(z~z') 



with the boundary conditions 



and 



cG prx (z,0) = d z ,G oS (z,z') 



z'=0 



G p<s (z,z')=0. 

Equations (|3.3[) can by solved as described in Ref. Ill], and one obtains 



G r(z,z') = 



exp (— \z — z'|) H — exp (— (z + z')) 



(3.3a) 



(3.3b) 



(3.3c) 



(3.4) 



where 
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Kj = ( T + p 2 + WA 2 



1/2 



(3.5) 



As in our previous work on RRN §,gjK)[|^{n| : this principal propagator may be viewed as being composed of a 
conducting and an insulating part: 



G p>x (z,z') = G™f{z,z')-G™{z,z') 



(3.6) 



where 



exp (-/^ |z - z'Q + ^ A ° exp (-k^ (z + z')) 
K X + c 



(3.7) 



and 



G™ (z,z') 



2Kg l 



exp (—Kg \z - z'|) H — exp (—Kg (z + z')) 

Kg + C 



^,0 



(3.8) 



B. Special transition 



Here we sketch our renormalization group improved perturbation calculation for the special transition. The ordinary 
transition is relegated to the next subsection. 

The constituting elements of our diagrammatic expansion are the vertex —g and the propagator G p ^ (z, z'). The 

decomposition of the principle propagator into G co S d (z, z') and G™ s (z, z') allows for a schematic decomposition of the 
principle Feynman diagrams into sums of conducting diagrams consisting of conducting and insulating propagators 
(see Fig. |^ and Fig. ^]). These conducting diagrams may be interpreted as being resistor networks themselves with 
conducting propagators corresponding to conductors and insulating propagators corresponding to open bonds. 
At the special transition G c °~ d (z, z') simplifies to the conducting Neumann propagator 



G c °f ,N (z, z') = ^— [exp (-k x \z - z'|) + exp (-k^ (z + z'))] 



(3.9) 



It is convenient to parametrize the conducting Neumann propagator in terms of modified Bessel functions of the 
second kind pfll: 



1 



P.A (4tt) 1/2 



x < exp 



ds 
TTTJ 



4s 



exp 



(r + p 2 



wX' 



exp 



(z + z'f 
4s 



(3.10) 



In this parametrization the parameters s correspond to resistances and the replica variables iXi to currents. The 
replica currents are conserved in each vertex and we may write Xi = Xi ^A, {^®1)> where X is an external current 

and denotes the set of independent loop currents. 



As we have learned in Sec. [ID, the surface coupling constant w$ is irrelevant. By setting it to zero, however, 



we would loose all information about its impact on the average resistance. To investigate this impact we work with 
insertions Here at the special transition we analyse correlation functions 



Gnm (jx,X||,^j ;r,i 
where the operator 



IN M 

D™ { ° s n v ( x *' n ^ ( x iU' 



H 



n - w 
Os-- 



/^SE( V ^ ( x ii'^)) 2 



(3.11) 



(3.12) 
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is inserted. Note by comparing Eqs. (2.25) and ( 3.1 2| ) that Os is associated with a coupling constant v$ — wg/w. 
To see how the insertion procedure works in detail we now consider diagram a (see Fig. |2j). Diagram a stands for 



a = — g w 



A(0 



E / ^ l) %% N (^)Gl% N (0,z'')G^ (z',z") 



where J is an abbreviation for (27r) ( d ^ j d d 1 p. The involved summation over the loop current reads 
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^TX^exp 



wP 



(a,aW) 



(3.13) 



(3.14) 



AC) 



We interpret P (X, XW) = - (s 2 + s 3 ) x^ 2 - 8l (a + aw) as the power of the diagram. To evaluate the summation 

we employ the saddle point method. Note that the saddle point equation is nothing more than the variation principle 
stated in Eq. fl2.7|). Thus, solving the saddle point equ ations is equivalent to determining the total resistance R({si}) 
of a diagram, and the saddle point evaluation of (3.14) yields 



(3.14) = — w exp — R (si, s 2 , S3) wX 2 



A(0_ 



91 



-X 



si + s 2 + s 3 



x exp 



(si + s 2 + s 3 )wX W5 



(3.15) 



Here we have switched back to continuous currents. The remaining integration is Gaussian and therefore straightfor- 
ward. In the limit D —>■ we obtain upon Taylor expansion of exp —R (si, s%, S3) wX 2 



(3.15) = -wX 2 



si 



si + s 2 + s 3 



O A 



(3.16) 



Note that is\X/ (s± + s 2 + S3) is the replica current flowing through the surface contact. Thus, we have an effective 
method to carry out the summations over the loop currents: essentialy we just need to determine the current flowing 
through the surface contact of the diagram. 

The ultraviolet divergent integrals occurring in the diagrammatic expansion can be regularized by the minimal 
subtraction method. We employ the following renormalization scheme: 



9S 



ryl/2 
Ifi — > Ifi = Z ' Ifi 

w — ► w = Z~ 1 Z w w 
+ <p s = Z 1 < Z 1 p s 



T — * T = Z 1 Z T T 

9 9 

c 



c = Zi 1 Z~ 1 Z c c , 



O s -» O s = Z 0s O s 



(3.17a) 
(3.17b) 
(3.17c) 
(3.17d) 



where G e — (4n)~ d / 2 r(l + e/2) with T denoting the Gamma function. Z, Z T , and Z u have been calculated to three- 
loop order by de Alcantara Bonfim et al ^2|. Diehl and Lam (l?]] computed Z\ and Z c to one loop order. Moreover, 
we gave a two- loop result for Z w in Hence, it remains to determine Zq s . Our calculation which is sketched in 
App. |A| yields 



los = 1 + 7- + O (u 2 ) 
6 e 



(3.18) 



The unrenormalized theory has to be independent of the length scale /i 1 introduced by the renormalization. In 
particular the unrenormalized correlation functions with the Os insertion have to be independent of /it, i.e., 



■S~Gn,m (^|x,X||, wX 2 | ;u,f,c, 



m -0. 



(3.19) 



Expressed in terms of renormalized quantities, Eq. (3.19) leads to the Gell-Mann-Low renormalization group equation 



N + M 



M 



.> -7 + y 7i + lo s 



G N ,M f{x,X||, wX 2 | 



; U, T, C, jJL 



, 



(3.20a) 







where we use the abbreviation 



d 







= n^r- + Ptt + TK ^r + W C^, — I- ck c — 

OfJ, ou 



d_ 

"dr 



d_ 

' dw 



d_ 



In Eq. ( 3.2C ) we have introduced the usual Wilson functions 

d 



7... («) = M 

C(u) =n 

K (u) = fl 



dfi 
d_ 

dfi 

d_ 

( \ d 
KS = fJ ''dfl 

0(u) = fi^ 



InZ 



In w = 7 — 7t, 



In t = 7 — 7r 



In c = 7 + 7i - 7 C 



= u (-e + 37 - 7„) 



(3.20b) 

(3.21a) 
(3.21b) 
(3.21c) 
(3.21d) 
(3.21e) 



where the bare quantities are kept fixed while taking the derivatives. 

The renormalization group equation can be solved in a standard manner by the method of characteristics. The 
characteristics read 



di 

At the infrared stable fixed point u* , determined by (3 (u*) — 0, we find 





with = fi , 


(3.22a) 




with u(l) = u , 


(3.22b) 


d 

I— lnr = k(w(0) 


with f(l) = t , 


(3.22c) 




with w(l) = w , 


(3.22d) 


d 

I— lnc = k s («(/)) 


with c(l) = c , 


(3.22e) 


InZ... = 7... (u(l)) 


with Z...(l) = 1 . 


(3.22f) 



G n ,m f {x, x|| , w\ 2 j;u,T, c, 



(Ar+M)j,/2+Mr ?1 /2+ 7 J. 



(3.23) 



where 7* = 7(1**), 7^ = 71 (u*), k* = = Kg (u*) and 7^ = 70 s (u*). In order to obtain the scaling 

behavior of the correlation functions, a dimensional analysis remains to be done. This dimensional analysis gives 



.fi 



(N+M)(d-2)/2+l 



G n ,m (jx,X||,wA 2 j ;it,T,e 

^|/ix, /iX|| , /i~ 2 uiA 2 j ;u* , (i~ 2 T, [i~ 1 c, /i 



xG 



N.M 



(3.24) 



Equation (3.23) in conjunction with Eq. (3.24) now results in 



G 



_ ,(Af+Af)(d-2+r,)/2+M J)1 /2+l+ 7 * e 



xG NM ({/x, lx h l- <t> / v w\ 2 \;u*,l- x ^T,l- 1 ^ s c, f j) 



Os 



(3.25) 



where 77 = 7*, 1/z/ = 2 — re*, 771 =7*, and l/vs — 1 — (to a re usual critical exponents for semi- infinite percolation, 
whose e expansion results may be inferred, e.g., from |22j|L7|. <f> = v (2 - £*) = 1 + e /42 + 4e 2 /3087 + O (e 3 ) is the 
bulk resistance exponent Q . 
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The scaling behavior of the correlation functions without i nserti on can be derived by similar means. The corre- 
sponding renormalization group equation is analogous to Eq. ( |3.20| ). Basically, just the 7e> 5 is missing. Solving this 
renormalization group equation gives in conjunction with dimensional analysis 



G NM ({x, X|| , W A 2 } ;u,t,c,h) = l(N+M)( d ~2 + v)/2 + M m /2 

x g nm ({fx, fx,, , r+i v w3? } ; u* , r rV us c , n) . 



By comparing Eqs. ( 3.25 ) and ( 3.26 ) we learn that the scaling dimension xq s of Os is 

XOs = 1 + l*Os > 

which tells us that the coupling vs scales in terms of the flow parameter / as 

V S {l) = V S l 1+1 °s . 

Taking into account that 

w(l) = wl C ~ 2 

we deduce that 

w s (l)=w s r' t ' s/u ■ 



(3.26) 
(3.27) 
(3.28) 
(3.29) 
(3.30) 



Here we introduced the surface resistance exponent <ps for the special transition. To one-loop order it is given by 

(3.31) 



*s = v{l-C-lh s ) = \-^+0{e 2 ) 



Now we can extract the scaling behavior of the average resistance between surface ports from the correlation 
functions with N = and M = 2. At the special transition, i.e., for r = c = 0, Eqs. ( |3.25| ) and (3.26) reduce to 



Go, 2 

and, respectively, 



■II ~ x ll 



G 



0,2 ( 



W \A =iW-2+^n)-W«-*)/"G 02 Ci 
) o s 'V 

, W A 2 ) = l^ d - 2 +^ Go 2 (i 



HI ~ x | 



HI ~ x | 



I _ II 



,r^X 2 ) 

l-t'^wX 2 ) , 



Os 



(3.32) 



(3.33) 



where 7711 = 77 + 771. Note that we have dropped some of the arguments for notational simplicity. Clearly, the two-point 
function with and the one without the insertion enter into the average resistance. Overall, we may write for the 
generating function 



G 



■II ~ x ll 



7A 2 ) =2( d - 2+f '»){Go,2 (l 

+ r (05-0)/^ Go 2 ^ 

We exploit the freedom to choose / = 



HI ~ x ll 



■II ~ x ll 



r*/»wx 2 



Os 



(3.34) 



HI - x || 



G 



HI ~ x ll 



,w\' 



HI - x || 



Moreover, we carry out a Taylor expansion which yields 

(d — 2+77,1 ) 



1 



•II ~ x ll 



WS 




(<j> s -(j>)/v 






x ll - x ll 




X ll " X || 


w 









wX 4 



wX 2 + 



(3.35) 



Here we have set all non universal constants equal to one. Upon differentiating Eq. (3.35) with respect to A 2 we 
obtain at A 2 = for the average resistance 

i>s—<f>)/v 



M «( x lh x ||) 







<t>/» r 

{ 1H 


w s 




j = w 


X|| - x|| 


w 


x ll - x ll 



(3.36) 



Note that (<j>s — 4>)/v < 0, i.e., the corresponding term is indeed only a correction to scaling that vanishes for 
X|| — xj| — > 00. 
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C. Ordinary transition 



At the ordinary tran sition the order parameter field satisfies the Dirichlet boundary condition <ps — 0. This can be 
deduced from Eq. (|3.4|), which shows that the Gaussian propagator vanishes for c — > oo if one of the points is located 
on the surface. Since correlation functions with insertions of the surface field are zero, they are not appropriate to 
gain information about fluctuation effects near the surface. A convenient method to investigate the scaling behavior 
of quantities that vanish for c — > oo is the 1/c expansion p3| , ^4f . For large c the propagator behaves as 



G Vl x(z,0) = c- 1 d z ,G° x (z,z') 



Olc 



z'=0 



(3.37) 



where G D - (z, z 1 ) = lim c ^.oo G n % (z, z') is the Dirichlet propagator. Similarly one can exand the entire correlation 
functions which gives 



G 



N,M 



(*H- x Il) 



c- M G 



N,M 



(Hi.*)) 



(3.38) 



where 



IN M 

G$ M ({x, X || , 9}) = Hm / J] <p (xi,$i) J] d n <p s (x, b ., 6, 



(3.39) 



3=1 



The subscript oo on the right hand side reminds us that the average is to be taken with c = oo. Apart from the 
additional term for N = 0, M = 2 the 1/c expansion amounts in replacing all surface fields tps by d n tps- The surface 
operator d n ips may be renormalized by the reparametrization 



dn^s -> [9n^] bare = (ZZ°Z) 1/2 d n <p s 



(3.40) 



A one-loop result for Zj 50 can be gleaned from pTy . 

We point out that Eq. ( ft.38| ) is correct in the given form for the bare, i.e., unrenormalized correlation fun ctions . 
In the renormalized theory nontrivial powers of 1/c will appear instead of 1/c. One may inferre from Eq. ( 2.15 ), 
however, that this peculiarity will not influence the average resistance. 

Here at the ordinary transition we insert the operator 



Of = | J d d - l x {] ]T (v s d nVs (x| h 



(3.41) 



into the one-loop diagrams contributing to G^ . The corresponding diagrams are depicted in Fig. |[ We use the 
renormalization 



Of -> Of = Z ~Of . 
Our one-loop calculation gives for the renormalization factor 

Details of the calculation can be found in appendix [b|. 

The renormalization group equation for the correlation functions G^ M with the Of insertion is given by 



(3.42) 



(3.43) 



N + M 



-7+Y7! 



= , 



(3.44a) 



where T>°° is an abbreviation for 



„™ 9 „ d d , d 

F Ofi ou or aw 



(3.44b) 



The Wilson functions 7f° and 7e>~ are defined as 
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7i°° (u) = li 



7a W (u) = /, ^ 



lnZo ? 



(3.45) 



_ (Af+A/)(d-2+J7)/2+Mr ; 5 x3 /2+3+7joo 



(3.46) 



Solving Eq. (3.44) provides us in conjunction with a dimensional analysis with the scaling behavior 

G N.M ({ X ' x ||> W ^ 2 } ; U ' T '^ 

x Gn.m (Ux, > r^^wX 2 } ; u*,r l ' v T, /it) 
where 77 J° = 7i° (it*) and 7^,00 = 70^ ("*)• For the correlation functions G^ ' M without insertion we obtain 

G% M ({x, x,| , w\ 2 1 ; u, r, M ) = i(N + M ){d -2 + ,)/2 + M vT ,2 

xG^^ M Uly.M\\,l^ lv w\ 2 Y,u*,l- l ' v T,^ . (3.47) 
By comparing Eqs. ( 3.46| ) and (3.47) we deduce that the coupling constant v^ 3 — wsJ (c 2 w) of Ol? scales as 

vf(l) ^vff^f . 

Hence we obtain for = ws /c 2 that 

Wf(l)=wfr*s' v , 

with the surface resistance exponent <jf^ for the ordinary transition reading 

4% = v(-i-c* 



1 19 

2 ~ 420 



(3.48) 
(3.49) 

(3.50) 



Now we can proceed in the same fashion as in Sec. [II B. We assemble the generating function for the average 
resistance from the two-point correlation functions (N = 0, M = 2) with and without insertion. Upon taking the 
derivative with respect to A 2 we finally arrive at 



M fl (x| h xf|) 



: ll _ x ll 



<f>/ v I 

1 + — 

w 



(3.51) 



D. Leading correction 

Now, as we have computed the corrections to scaling due the resistors located on the surface it is legitimate to ask: 
How important are these corrections compared to other corrections? In particular we should compare them to the 
leading correction which is governed by the so called Wegner exponent. This leading correction emerges when the 
renormalized coupling u is not exactly equal u* since the renormalization flow has not arrived at its fixed point yet. 
Such a case occurs typically when there is a finite momentum cutoff reminiscent of a no nvanishing lattice spacing. 

In order to determine the leading correction we revisit the characteristic Eq. (3.22b). Upon expansion for small 
deviations u — u* from u* we obtain 



I— = us [u — u*\ + O (ju — u*] 2 ^ with u(l) = 

where u> = /3' (u*). This differential equation is readily solved with the result 

u(l) = u* + [u - it*] l u . 



(3.52) 



(3.53) 



lu is referred to as the Wegner exponent. It can be calculated without much effort to third order in e upon using the 
three-loop result for (3 (u) obtained by de Alcantara Bonfim et al. |2^] . Here we are working only to first order in e to 
which the Wegner exponent is given by 
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O e 2 



(3.54) 



To see the impact of the deviation from u* on the correlation functions, we revisit the special transition. The 
analysis can be adapted, however, in an obvious fashion to the ordinary transition. From the renormalization group 
equation and Eq. (3.53) we deduce that 



Go 



■II 



M 2 ) =i {d - 2+ ^G , 2 (i 

x j 1 + [u - u*] l u F ( I 



wX 2 + O [u - u 



*i2 



(3.55) 



where F is a scaling function. Once more we choose I 

(d-2+J7|,) 



G 



0,2 ( 



,wX- 



and carry out a Taylor expansion, 

4,/u 



wX 2 + 



x \l + [u-u* 



1 



wX 2 + 



(3.56) 



where we have set, as always, all non universal expansion coefficients equal to one. By taking the derivative with 
respect to X 2 we obtain at A 2 = 



M 



Jl(x||,xj|) 

















^ = w 


x ll - x ll 






- [it — u*] 


x,| - x|| 





(3.57) 



Now we can compare the different corrections to Mr considered in this paper. To first order in e the correction due 
to the deviation from u* falls off algebraically for increasing port separation with an exponent — lu = — e. The surface 
correction at the special transition vanishes much faster with an exponent (<^s — <fi)/v = —1 + e/21. Among the three 
corrections, the surface correction at the ordinary transition drops off fastest with (</>|? — 4>)/v = — 3 + 23e/105. 



IV. CONCLUSIONS 



In order to study the effects of surfaces on resistor percolation we have considered a semi-infinte RRN. We have 
presented a field theoretic Hamiltonian in which the coupling constant corresponding to the surface conductances 
turned out to be irrelevant. We have calculated the corrections to scaling due to this irrelevant coupling for the 
special and the ordinary transition to one-loop order. 

In this paper we did not consider the surface and the extraordinary transition. We left the surface transition aside, 
because it is basically equivalent to the percolation transition of a translationally invariant d — 1 dimensional RRN. 
Thus, the behavior of Mr at the surface transition can be inferred from Ref. Q . We did not drill into the extraordinary 
transition because of severe technical complications. These are rooted in the fact that the order parameter profile is 
not flat on neither side of the line of the extraordinary transition. 

That the surface coupling ws is irrelevant seems intuitively plausible. Suppose that the resistor network is at 
the special transition where the surface and the bulk percolate simultaneously. Assume that we conduct a series 
of consecutive measurements in which we apply the external current / between two surface ports x» and x'« with 



increasing distance 



As we increase 



more and more paths of connected bulk resistors will add in 



parallel to the connected paths of surface resistors. Hence, the influence of the surface bonds becomes negligible for 
large 



pcrc 



This situation is even more pronounced at the ordinary transition, where the percolation probability 

is lower at the surface than in the bulk. Thus, it is plausible that the surface correction to Mr vanishes faster 

for increasing x\\ — x'« at the ordinary transition than at the special transition. 

Both surface corrections turn out to be small compared to the leading correction governed by the Wegner exponent 
for percolation. This means that the surface has weak effects on the average resistance compared to those of a finite 
lattice spacing. As long as one is interested only in the leading behavior and mayor corrections to it, one may neglect 
the surface effects safely. 

To our knowledge there are no numerical simulations available to date which could be used to test our predictions 
for the surface resistance exponents. The reason is probably that the typical simulations reported in the literature 
use a so called bus-bar geometry, in which the resistor network is placed between superconducting plates which short 
entire surfaces. We hope that this paper triggers simulations providing numerical estimates for </>s and 0§? . 
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APPENDIX A: EVALUATION OF DIAGRAMS FOR THE SPECIAL TRANSITION 



In this appendix we sketch the computation of Zq s for the special t ransi tion. We start with diagram a and revisit 
Eq. (3.13). Upon inserting the result for the current summation, Eq. ( 3.16| ), we obtain 



a = -g 2 w\ 2 



1 



dsid,S2ds3 



Sl 



x exp 



(47r) 1 Jp Jo VS1S2S3 \si + s 2 + s 

,2 



exp [- (si + s 2 + S3) (t + p 2 )] 



4s 2 4s a 



exp 



(z-z'Y 
4si 



exp 



(z + z'Y 
4s 1 



(Al) 



where we have dropped all other terms since we are interested here only in the part of a proportional to w. The 
momentum integration is straightforward and yields 



a = -g 2 w\ 2 



1 



dsid,S2ds3 



si 



(47r) 3/2 Jo V S 1 S 2S3 V s i+ s 2 + s 3 / V 47r ( s i + s 2 + s 3 ) 



1 



x exp [- (si + s 2 + s 3 ) t] exp 



4s 2 4s a 



exp 



(*-*y 

4*i 



exp 



(d-2)/2 



(z + zf 



4sx 



(A2) 



At this stage it is useful to apply the Laplace transformation. The Laplace transformed of a reads 



/>oo poo 

C (a) = / / dz' exp [— uz — vz'\ a 

Jo A) 



= -9 2 w\ 2 



(si + s 2 + S3 



(47T) d/2 JO 
o2 

1 ,2+d/2 + 



dsids 2 ds 3 exp [— (si + s 2 + S3) r] 



(A3) 



where we dropped higher order terms in u and v since they are convergent. The integrations in Eq. (A3) are simplified 
by the change of variables si — > tx, s 2 — > ty, and S3 — > t (1 — x — y): 



C (a) = -g 2 w\ 2 



(4tt 



d/2 



r-1 fX—X pOQ 

dx dy dt x 2 t 2 ~ d/2 exp [-tr] 

Jo Jo Jo 



The result 



C (a) = -g 2 w\ 2 
can be expanded for small e = 6 — d as 



1 V/2-3 r U_d 

(4^ 2 6 ^ 2 



£ ( a ) = -g 2 w \ 2 r-^ 2 



3e 



Upon transforming back to real space we obtain 



a=-. 9 2 ^A 2 ^T-/ 2 <5 + (z ' )(5+ 



(A4) 



(A5) 



(A6) 



(A7) 



with <5_|_ (z) denoting the distribution defined by 
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6+ (z)g(z)=g(0) 



Alternatively, a can be computed with help of the parameter sum 



E (n, to, k) = lim V A (/) 
The evaluation of this sum yielding 



1 



K X+\m K X v^a+aco + k a 



(A8) 



(A9) 



E { n, to, *) = A 2 1 2- "(" + 2) + (n + l)fc + ^(fc-l) 
v ; K 'i +m+fe (n + to + k) (n + to + k + 2) 



+ A 



(A10) 



is outlined in appendix In this approach we do not parametrize the Neumann propagator. Instead, we substitute 
Eq. ( |3.9[ ) directly into Eq. ( |3 .13 ). Laplace transformation then gives 



C (a) = -g 2 w 



+ 



AW 

1 



|_ A 

1 



2k^ + u + u \ K X+Xw +^x + u K U $ W + K Z + V 



(All) 



Taylor expansion of the right hand side of Eq. (All) leads to 



£(a)=- g ^E/ X( ° % 1 2 



K A+A<0 + K A 



+ 0(u,v) 



(A12) 



In terms of the parameter sum, Eq. ( A12 ) reads 



1 



£ (a) = -g 2 w- / E (1, 3, 1) + E (1, 2, 2) + O (n + m + k > 5) 



The terms with n + to + fc > 5 are convergent. We keep the divergent contributions and obtain 



(A13) 



(A14) 



The remaining momentum integration is straightforward. In real space one retrieves the result stated in Eq. (A7) 
The second diagram of the right-hand side of FIG. 0, b, is easier to compute. We mention only the result 



e 



IH+ (z')6 + [z") 



(A15) 



APPENDIX B: EVALUATION OF DIAGRAMS FOR THE ORDINARY TRANSITION 



We start with diagram c displayed in FIG. ||. It stands for 



AW 



2 wJ2 / A ( ' )2 



a;=0 



(Bl) 



where the conducting Dirichlet propagator reads 



1 



G d1 ' ( Z ' Z ') = 9~ [° XP \ Z - Z '\) - CX P (- K X ( Z + Z '))\ 

P. LK^ 



(B2) 
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Upon insertion of Eq. (B2) into Eq. (Bl) one finds for the Laplace transformed of Eq. (Bl) 

1 ( 1 1 



\+\( 



(B3) 



Next we carry out a Taylor expansion in terms of 1/kj and 1/ (^ X(!)+x + n^j • We keep only those terms proportional 

to uv, since these are giving the leading behavior in the limit c — > oo. They are in real space proportional to 
S' + (z 1 ) 6' + (z"), where the distribution S' + (z) is defined by 



5>(z)g(z)=g'(0) . 



We obtain in terms of the parameter sum E (n, m, k): 

C (c) = —g 2 wuv 
Transforming back to real space yields 



iE(l,3,l) + ^S(l,2,2)-lE(l I 0,4) + O(n + m + fc>5) 



~ g2wP We T ~ t/2S '+ {Z ' ] 6>+ {Z " ] 



For the second diagram on the right-hand side of FIG. |^ we obtain by similar means 



-q 2 wX 2 r 

y 6e 



5G £ 



^ 2 S'(z')5'(z") 



(B4) 



(B5) 



(B6) 



(B7) 



APPENDIX C: EVALUATION OF THE PARAMETER SUM 

Here we outline the evaluation of the parameter sum introduced in appendix [A|. We apply the inverse Mellin 
transformation [G5| 



da r (a) T(k — a) _ (T _ 

r(fc) 



for < a < k, to the right-hand side of Eq. (A9). Using Schwinger parametrization we obtain 



E (n, m, k) — lim 

D->0 



da T (a) T{k-a) 



dsi / ds2s\ 



(n+a)/2-l {m+k-a)/2-l 



ioo 2 7rir(^)r(2^^)r(fc) J L J Q 

exp[- (si +s 2 ) (r + p)] \-—q— ) X] exp - s i w - s 2 w\ {l)2 



AW 



(CI) 



(C2) 



Upon completion of the square the sum is easily carried out in the limit D — > 0. Moreover, we expand for small A 2 . 
We obtain 



E (n, to, k) = A z 



da T (a) T(k-a) 



OL — too 
2 



27 rir(^)r(2^^)r(fc)7 Jo 



oo p oo 



dsi / ds2s\ 



(n+a)/2-l {m+k-cr)/2-l 



(si + S 2 ) 



^—2 exp [- (si + s 2 ) (t + p)] , 



(C3) 



up to terms of higher order in A 2 . A change of variables si — > tx, S2 — > i(l — x) renders the integration over the 
Schwinger parameters straightforward. One gets 
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£ (n, m, fc) = A z 



1 



da T (cr) T(k-a) (n + a) (n + a + 2) 



(n+m+l) I 0_„- 



r (fc) (n + m + k)(n + m + k + 2) 
The remaining integration can be done by exploiting the identity 



da T (a) T (fc - cr) v 
-a 



ioo 2m T(fc) V dl 

Finally, one arrives at the result stated in Eq. ( A1C| ) 



d_V r a - io ° da T{a)T(k-a) t _ 



2ni 



T(fc) 



(C4) 



(C5) 



t=i 
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FIG. 1. Schematic phase diagram for semi-infinite percolation. The horizontal axis corresponds to — r ~ p — p c . The vertical 
axis corresponds to c, whose negative is a measure of the surface enhancement. The lines labeled ordinary, extraordinary and 
surface indicate continuous phase transitions that have been given these names by Lubensky and Rubin Jl!^]. The lines meet 
at a tricritical point that represents the special transition. 




FIG. 2. Ultraviolet divergent Feynman diagrams contributing to the renormalization of Os at the special transition. The 
bold lines represent principal propagators G p jj (z,z'). The light lines stand for conducting propag ators G c °l A {z,z') and the 

dashed lines represent insulating propagators G™ s (z, z'). Insertions of the surface operator Os are symbolized by the hatched 
" surface" . 




FIG. 3. Ultraviolet divergent Feynman diagrams contributing to the renormalization of 0%° at the ordinary transition. The 
meaning of the symbols is the same as in FIG. [| except that the hatched "surface" stands now, in conjunction with the bars, 
for an insertion of the operator 0%° ■ 
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